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Abstract: We study the functional dependence of the spin- weighted angular mo- 
ments of the two-point correlation function of the three dimensional cosmic shear on 
the expansion history of the universe. We first express the redshift dependent total 
equation of state parameter in terms of the growing mode of the gauge invariant 
metric perturbation in the conformal-Newtonian gauge for the case of adiabatic per- 
turbations. We then express the redshift dependent angular moments of the shear 
two-point correlation function as an integral in terms of the metric perturbation. We 
present the final explicit expression for the case of a Harrison-Zeldovich spectrum of 
primordial perturbations. Our analysis is restricted to the linear regime. We use our 
results to make a preliminary study of the required sensitivity that will allow cosmic 
shear observations to add significant information about the expansion history of the 
universe. 
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1. Introduction 



Determining the expansion history of the universe is one of the central problems 
in cosmology, and the goal of many observation programs: distant supernovae [|l], 
1^, 1^, the large scale structure in the universe 0, ||, and the cosmic microwave 
background 0. It has become clear in recent years that kinematic distance probes 
in the homogeneous and isotropic universe, such as luminosity distance, angular 
distance etc., have a limited resolving power for determining the expansion history 
of the universe. They rely on light emitted from distant sources, and hence they 
measure an integral over the expansion history [^, ^, |TD[ . It is possible to improve the 
determination of the expansion history by adding prior assumptions on the evolution, 
or by combining several kinds of observations ||lT], |l2l, or by focusing on better 
determined quantities O, 111, llH]. It has also been suggested that the measured 



integrals can be differentiated |jT6[. Reviews of the subject can be found in IjT^ and 



m 



Cosmological perturbations provide, through their dependence on the homoge- 
neous and isotropic background, an independent tool to probe the expansion history 
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of the universe |19, 20, 21 1. To be able to use the perturbations to determine the 
expansion history of the universe it is necessary to express the total equation of state 
parameter Wtot in terms of the perturbations in a way that does not depend on the 
specific functional dependence of Wtot on time. In the following we refer to such an 
expression as "model-independent". Further, it is necessary to identify observable 
quantities from which one can determine in a reliable way the perturbations as a 
function of time, or equivalently, of redshift. Obviously, we have to look for ob- 
servables that can be measured precisely, however, in addition they also have to be 
evaluated precisely, otherwise the theoretical errors form the imprecise calculation 
will dominate the final error budget and will limit their resolving power as probes of 
the expansion history of the universe. 

The three dimensional cosmic shear ||22|, ^ |2^ seems to be a promising observ- 
able that can be measured by weak lensing observations. Currently, the observations 
are mostly two dimensional [p5| , |2^ , p7| with preliminary studies of 3D analysis |]28|, ^ 
and future programs expected to have three dimensional capabilities PDI The 



recent advances in measuring the cosmic shear and the expected future improvements 
have attracted a growing interest in the potential of weak lensing measurements for 
determining the expansion history of the universe, either by itself ^2 



I 34, 35, 36 



or in combination with other measurements ^ Cosmic shear is a measure 
of the shape distortion of distant objects as the light emitted from them propagates 
through the perturbed universe. To measure the three dimensional cosmic shear the 
redshift of the distorted objects (usually galaxies) needs to be measured in addition 
to the distortion pattern. Several recent reviews of weak lensing give a comprehensive 
and exhaustive description of the current state of the field ^ . 

If the perturbations are weak, the cosmic shear depends linearly on them. In this 
linear regime of the metric perturbations, both observations and theoretical calcula- 
tions can be done reliably without any obvious obstruction. Thus, it seems possible 
to reduce the intrinsic theoretical and experimental errors down to the percent level, 
as required for accurate determination of the expansion history of the universe. In 
gauge invariant perturbation theory, the metric perturbation $ can in principle be 
small even if the density perturbations Sp/p are not, thus the use of metric pertur- 
bations allows us to extend the linear regime to smaller scales. 

The standard approach of exploring the relationship between the cosmic shear 
and the expansion history is to use numerical methods. Numerical comparison to 
cosmological models is made and used to estimate the prospects of constraining 
cosmological parameters from shear measurements. This method gives precise results 
for specific models and for specific parametrizations. In this approach the main 
important factors that need to be determined are related to the accuracy of the 



shear measurements and to the estimators for the power spectrum |44 



We wish to take a different approach. We pose a different question which we 
believe is quite significant and moreover can be answered in a definite way. The 
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question that we wish to pose and answer is: What is the theoretical hmit on model- 
independent information about the expansion history of the universe that can be 
obtained from cosmic shear measurements with a given accuracy? We wish to know, 
in theory, whether the 3D cosmic shear (or its angular power spectrum) is more 
sensitive to changes in the total equation of state parameter than other observables. 
Our results allow us to estimate the accuracy goal needed for shear measurements so 
they can improve on other accurate tests such as luminosity distance measurements 
or CMB measurements. Our approach is mostly relevant when trying to estimate 
the prospects of future lensing surveys for constraining the evolution of the universe 
in a model independent way. To answer the question that we have just posed we 
need to determine the functional dependence of the cosmic shear on the expansion 
history of the universe. In the context of this paper this is equivalent to determining 
the sensitivity of the angular spectrum of the 3D cosmic shear to changes in the total 
equation of state parameter. 

In this paper we find a model independent analytical solution for the growth of 
the metric perturbations (Section |^) and show that the spatial and time dependent 
parts are separable under certain conditions. We then use this solution to show the 
functional dependence of the three dimensional cosmic shear on the DE equation of 
state parameter (Section ^) in the linear regime. We then explain how to extract 
Wtot from the red-shift evolution of the shear angular multipoles (Section ^). Using 
our results we make a preliminary analysis of the theoretical sensitivity to changes 
in the evolution of dark energy and estimate the influence of such changes on shape 
and strength of the shear angular spectrum. Section ^ contains our conclusions. In 
appendix ^ we present a detailed derivation of the functional dependence of the total 
equation of state parameter on the metric perturbation. 

2. Expressing $ in terms of Wtot 
2.1 The perturbation equations 

The line element of the perturbed universe in the conformal Newtonian (longitudinal) 
gauge is 

ds^ = a^{r]) [-(1 + $ {r], x)) drf + (1 - $ [t], x)) [dw^ + fl{w)dn'^)] . (2.1) 

7] is the conformal time, dVL is the two sphere differential element, and fx depends 
on the spacial curvature 

( K-^l^ sin{K^'^w) ,K>0 
fK{w) = 1 w ,K = 0. (2.2) 

[ (-ir)-i/2sinh((-ir)i/2«;) ,K<0 

We choose to ignore shear perturbations and assume that the cosmic fluid is a perfect 
fluid. Imperfect fluid perturbations' influence is in general considered to be small 
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(see for example and could be evaluated in subsequent investigations. We follow 
the standard derivations that are reviewed in |^ to obtain the equations of motion 
for $. They are derived from the perturbed Einstein's equations, 

- - 3(7^^ - K)^ = ATT\K\Ga'^6p, (2.3) 
$" + 37^$' + {2n' + {n'^ -K))(!> = A7r\K\Ga^6p. (2.4) 

where Ti = a' /a, and prime denotes a derivative with respect to rj. The limit K —* 
of these equations has to be taken carefully keeping Ka^ fixed, 

As mentioned above, the linear equations are valid when $ is small |$| ^ 1, 
however this does not require in general that 6p/p is small. From eq. (|2.3|) we can 
estimate that for small wavelengths 5p/p is larger than $ by a factor of the order of 
the square of the ratio of the size of the horizon to the wavelength, Sp/p ~ (g/7i)^$. 

Substituting 6p = {^)s^P + {^)p^S = c^Sp + t6S, where the total speed of 
sound of the perturbations is c| = ^ and 6S is the total entropy perturbation, leads 
to a single second order equation for 

+ 3n{l + 4)$' - 4V^$ + {2n' + (l + 34)(^^ - K))<!> = A7i\K\Ga^T6S. (2.5) 

In the rest of the paper we will only consider adiabatic perturbations for which 
6S = 0. Then, 

$" + 3n{i + 4)^' - 4v^$ + [2n' + (i + 34){n'^ -k)]^ = o. (2.6) 

Here we can take the K ^ limit in a straightforward manner. 

Our derivation is fully relativistic. In doing so we can put initial conditions on 
$ outside the horizon and follow its evolution, and hence use directly the early-time 
information about the spectrum of metric perturbations from the CMB or the linear 
matter power spectrum rather than using the late time processed matter power spec- 
trum. The numerical difference between the relativistic and non-relativistic analysis 
at small redshift z for small wavelength perturbations is not expected to be large. 
Again, we can roughly estimate the difference from eq. (|2.3|) to be of the order of the 
square of the ratio of the size of the horizon to the wavelength ~ [q/Ti.)^. 

In the perturbed Einstein equations ( |2.3| - |2.4| ) $ is related to the total density 



perturbation 6 p. In a model of a universe containing dark energy and cold matter the 
solution for $ depends on the two component background and on both perturbations, 
the dark energy perturbation and the cold matter perturbation. Equation ( p.6|) is 
always correct when the total speed of sound c^^ = is used. If the various 

components are weakly coupled as expected for matter and DE, then 
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and thus the total speed of sound is {csf'Sptot = J2i {(^Sif^Pi- For the two fluid 
model 

{csf^ptot = {csmfSpm + {csdeY^Pde- (2.8) 

Here the subscript m denotes matter quantities and the subscript DE denotes DE 
quantities. The matter speed of sound vanishes {csm)"^ = 0, thus, for adiabatic 
perturbations we get 



$" + 31 + ^-f^icsnEY n^' - -f^icsDEYV'^ 



+ 



2n'+[i + ?>^-p^{csnEf) ire-K) 

Optot J 



$ = 0. (2.9) 



2.2 The solution of the perturbation equations 

The general solution of eq. (|2.6|) is conveniently obtained by a standard change of 
variables to m, 

^ = AnG^/^Tpu. (2.10) 



The equation for u is 



with 



nil 

u" -clV\- ju = 0, (2.11) 



Here we assume spatial flatness. The expression and solutions for curved space can 
be found in |^6|]. Since the in the universe space curvature is known to be quite small 



we expect to be able to treat it as a perturbation in subsequent analysis. 

We can express 9 in terms of the time-dependent total equation of state param- 
eter 

wtot = -, (2.13) 
P 

Finding the full solution of the perturbation equations requires solving the two 
fluid equations. However, we are interested in the growing solution at rather late 
times (say 2; < 4), when substantial deviations from matter domination start to 
build up and for wavelengths that are smaller than the horizon. At those late times 
the perturbations in the cold matter will be the dominant perturbations and we 
will be able to safely ignore the DE perturbation while taking into account the 
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changes in the background evolution due to DE (see [47, 48] for a recent discussion). 
As we explain below, if the perturbations are adiabatic then the initial amplitudes 
of the different types of perturbations are about equal at horizon entry. During 
matter domination (when Wtot ~ 0) the well known "growing mode" solution of 
the perturbation equation is constant. DE perturbations, on the other hand, decay 
during matter domination. 

To understand our argument more precisely, let us consider the following situ- 
ation. Let us assume, for the moment, that the DE perturbation dominates. The 
general solution assuming that the DE perturbation is the dominant one in Fourier 
space, for a mode with wave vector q, is 



u = VviCiJ^dV^v) + C2Y^/2{V^QV))- (2.15) 

For large rj (late times) the Bessel functions decay as s/fj. Consequently, the 
solution for u approaches a constant at late times. From eq.( |2.1(JD , since during 
matter domination p ~ and a{ri) ~ yf , it follows that the solution for $ decays 
as T)'^ 

$ ~ t]-^u{t]) ~ 7]-^. (2.16) 

We see that if the DE and matter perturbations start off with equal amplitudes, the 
DE perturbation will decay through the era of matter domination with respect to the 
matter perturbation by a factor of ri~^ = a^/^ = (1 + z)"^^"^. For example, at z = 1 
a DE perturbation that entered the horizon at z = 100 will be smaller by a factor 
of about 10~^ than a matter perturbation that entered the horizon at the same time 
with the same amplitude. 

Now, let us focus on the matter perturbations. We can solve the equation for 
the matter perturbations in (|2.9|) with an arbitrary background equation of state 



Wtot- As the relative part of the dark energy goes to zero the value of the total 
speed of sound c| goes to zero. This is equivalent to the equation for a single fluid 
with a vanishing or negligibly small c|. The exact condition on eq.( |2.11| ) that we 
will assume in Fourier space, for a mode with wave vector q is 

{qvfcl « 1. (2.17) 

Notice that we restrict ourselves to a positive speed of sound for the dark component. 
Although a negative speed of sound isn't prohibited these solutions are unstable and 



restricted in the DE's equation of state parameter space. The solution of eq. ( |2.11 
is given by 

u{x, ri) = Ci{x)e{ri) + C2{x)e{'n) [ d7j-^. (2.18 



The first term is the smaller decaying solution and the second is the larger term 
which is usually referred to as the "growing solution" even though it is sometimes 
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constant, or decays slower than the first term. We are interested only in the growing 
solution, because it dominates the solution at late times. 

Using eq.( |2.14| ) we can express the growing solution of eq.( |2.18| ). 



$+(r/,a-) =47rGC2(x)^ / d?ja\?j){l + w,,,{v))- (2.19) 

a J 

From eq. (|2.19|) we can see that the contribution to $+ from an era when Wtot = 
— 1 vanishes. This is to be expected since Wtot approaches a constant value of — 1 
only if the DE is a cosmological constant and it does so at very late times when 
the cosmological constant completely dominates the matter. The solution x) 
factorizes into a time dependent and spatial part 

$+(//, f) = C(x)$T(r/) (2.20) 
C{x) = A7fGC2{x) (2.21) 

^Av) = — [ d?ja\?f){l + wUv)) (2.22) 



The time dependent part $t obeys the following differential equation in redshift 
space 

<^t{z) - ^^^^""l , d,<^T{z) = (2.23) 
Solving for lUtot in terms of $ we find 

The details of the derivation of eqs. (|2.23|) and ( p. 241) can be found in Appendix |^. 

The relation between the density perturbation and the metric perturbation is 
determined (for a spatially fiat universe) by the following equation 

^ = -1^ rv2$ - 3W - 3n^^] . (2.25) 

For the ACDM model, if one uses the solution for ^t{z) it is possible to show that 
it is identical to the known solution for the linear density perturbation growth factor 
which was obtained using different methods p9| , 

oo g 

$r"(^) = + / ^^'(1 + (^) • (2-26) 

z 

The details of the derivation are given in appendix 0. Our solution for ^t{z) there- 
fore amounts to a generalization of the known solutions for the linear growth factor 
to the case of an arbitrary equation of state of the DE. 
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So far we have not discussed the space dependent factor C{x), which we do 
now. The primordial spectrum of C{x) is an input for our analysis. It is usually 
assumed to be a power-law spectrum, and that the perturbations are isotropic and 
homogeneous. The primordial spectrum of $ can be parametrized as 

P^{k) = A{kr]or-' (2.27) 

The parameter n is the spectral index and A is the spectral amplitude. Both were 
measured most recently by the WMAP experiment 0. In particular, the spectral 
index is approximately n = .95 corresponding to an approximately flat spectrum. 
The range over which the spectrum is flat {n is approximately 1) is limited because 
causal processes inside the horizon suppress the perturbations [^ . 

The solution for $ above is valid only after matter domination. So the primordial 
spectrum has to be evolved into the "initial" spectrum at the beginning of matter 
domination. We shall use for this purpose the standard practice of including a 
transfer function Tk{ri), 

^k,r])=TkivMk). (2.28) 

The transfer function is normalized such that Tk{i]o) = 1 for k —>■ 0. The two-point 
correlation function is then given by 

($(A;,r7)$ * {k',r]')) = Tk{v)Tk'W)k-'P'fik){2nf6%k - k') 

= n{7])Tkiv'){^7r)'{kVor-'A6'{k-k'). (2.29) 



To obtain the value of the initial spectrum we must input into eq. ( |2.29| ) the value of 



1] at the beginning of matter domination- rjin. The initial spectrum is then given by 



ZTT \ / i 



{C{q),C*{q'))=T,\v.r.){- A6{q-q'). (2.30) 

The value of Tq{j]ir^ can be evaluated using several approximations or calculated 
numerically using CMBFAST/CAMB etc.. Since during matter domination the per- 
turbations are frozen ($ is constant) the exact value of is not of particular im- 
portance. 

3. Expressing 7 in terms of $ 

Although it is common to write the shear as a function of the lensing potential we 
choose to leave it as a function of the metric perturbation using the solution presented 
in SecH The expression for the shear is 



7 = 7i + «72 = / — 



+ 2i 
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For a single source at distance w the shear is given in eg . (13 .11 ). We would like even- 
tually to find the angular moments of the shear-shear two-point correlation func- 
tion. Because the shear is not a scalar we have to use the spin-weight spherical 



harmonics formalism as in El]. The relevant properties of the s- weighted spheri- 



cal harmonics syi,m can be found in [pTl p3. The spin- weight operator 9 operates 



on the s-weight spherical harmonic and gives an s + 1-weight spherical harmonic 
9 s^/,m = [(/ — + s + s+i^,m- Exprcsscd in terms of the spin- weight 

operator the shear is given by 

^{w,e,^)= fdw' ;f] m<!>iw',e,^). (3.2) 

J fK[W)fK{w') 


For a spatially fiat universe fxiw) = w as can be seen from eg. (|2.2| ). Recall in addi- 
tion that $ can be factored into a space dependent factor C{x) and time dependent 
factor $T, defined in eg. ( |A.8| ). Combining these facts we arrive at the final expression 
for the shear, 



,w — w' 



7(w, 9, ^) = dw' <l>T(w')m C(w', e, (3.3) 

WW 







Because this is a spin-weight 2 object it can be decomposed into an even and odd 
parts, these correspond to the well known E and B modes. Using the fact that the 
B modes of the shear field vanish and to keep things simple we will compute the 
correlation for the full expression. A detailed explanation of the decomposition and 
its properties can be found in [^. Now, let us compute the two-point correlation 
function, {■j{wi,9i,(pi),'j*{w2,92,(p2))- We can use the standard Fourier expansion 
and the assumption that (C(gi), C*(g2)) = f{<l)S{(li — ^2) to obtain 

(7(^1, 6*1,(^1), 7* (^2, 6*2, V52)) = 
Wl W2 

= fdu, fdu2^^^^^^^^^^<^Tiui)^Tiu2)QiQi32^2 fd\ f {q)e''^<^^~^^'\i3A) 

J J Wl Ml W2 U2 J 



We now expand the exponential e^^'^^'^~^^^ in ordinary spherical harmonics and per- 
form the integration on the unit sphere in g-space, 

{liWl,9l,(pi),-f*{w2,62,(p2)) = 

Wl W2 

= {4nf fdu, ^M, ^^~^^^^~^^ $H^i)$H^,) fdqq'fiq) 

J J WiUi W2 U2 J 



^3i{<iui)ji{qu2)^{di Yim{9i, <^i)B2B2 yL{02, V^2). (3.5) 





X 

Lm 



The two spin-weight operators act on the ordinary (spin-weight zero) spherical har- 
monics in the expansion and give spin-weight s = 2 spherical harmonics, and similarly 



- 9 - 



the conjugate spin- weight operators give spin-weight s = 



—2 spherical harmonics, 



Wl W2 



(4vr)2 [du, [du2^^^^^^^^^^^T{ur)^T{u2) [dqq'f{q) 

J J WiUi W2 U2 J 



X 





^ (1-2)! -^''^^^^^-^''^^^^^ 2>^m(^l, Vl) -2Yi*^{e2, <^2)- (3.6) 



We may perform the summation over m using the summation rule for spin-weight 
spherical harmonics 



^ siY,:^{ei,y^i) S2YU()2,V2) = y ^ .,ll,_.,(Aa)e-^^ (3.7) 

m 

The angels a,P and 6 are the rotation angels from {9i,{pi) to {02,(p2)- The two- 
point correlation function should only depend on (3, being the angle between the two 
directions. Hence we may choose the polar axis of (f such that it is aligned with the 
two points and set a = 5 = 0. In this case, our final expression for the two-point 
correlation function is 

iliwi, 6*1, V9i), 7*(u;2, 02, V2)) = 

= {Airf fdu, [du2^^^^^^^^^^<^Tiu^)Mu2) fdqq'fiq) 

J J WiUi W2 U2 J 





X 



E V^f3||-^'(^"i)^'(^"2) 2Y1MO). (3.8) 



The shear spin-weight 2 angular power spectrum 

The shear can be expanded in the s = 2 spin-weight spherical harmonics 7 = 
J2im 20-im 2Yim- From the definition of 7 in eq.( ^.2| ) it follows that 7* is proportional 
to -2^m- The conjugation relation of spin-weight spherical harmonics implies that 
2aim = -20'im- From the isotropy and homogeneity of the shear-shear two-point cor- 
relation function we know that it must be a function of |x — x'| only, so the two-point 
function of the coefficients 2fl«m can only depend on /, ( 2a«m, 2(^1 m') — ^u'^mm' 2C1. 
Consequently, we may express the shear-shear two-point function in terms of the 

angular spin-weight two coefficients (7,7*) = Xlz \J 2C1 2^z,-2- The summation 
over m was performed using eq.( |3.7|) . By using the orthogonality relationship of 
the spin- weight spherical harmonics / dVL gYim sYvm' — ^ii'^mm' we can extract the 
angular coefficients 2C1 from eq. (p.8|) 

Wl W2 

2Q = (4vr)^^|^ [du, /rftX2 '^'~"''^'~"' $T(^i)$T(^2) fagV(g)j/(g^i)j/(g^2i3.9) 
{I -2)1 J J wiui W2U2 J 
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Figure 1: The shear spin- weight 2 angular power spectrum for different red shifts. The 
cosmological model is ACDM with {^Ide = 0.7, ^Ik = 0,n = 1,wde = The four hues 
correspond to redshifts z = 0.4, 1, 2, 3 from the bottom up, respectively. 



We can also recall now that for a flat spectrum (n = 1) the value of f{q) is given 
in equation (|2.30| ) so that 

fiq)=A (yyr/(r/,„). (3.10) 

The /-dependence of the multipoles of the cosmic shear angular power spectrum 
for the case n = 1 is shown in Fig. |I] for three values of redshift. The redshift 
dependence of the multipoles of the cosmic shear angular power spectrum is shown 
in Fig. ^ for three values of /. 

The redshift dependence of the multipoles of the cosmic shear angular power 
spectrum is shown in Fig. ^ for three values of /. 

The expression for the shear spin-weight 2 angular power spectrum in eq. |3.9| has 
a highly oscillatory integrand due to the factor of the spherical Bessel functions. For 
the case of Wi = W2 we know that for every q there is a cutoff distance such that 
qu ~ / where the integrand practically vanishes. Since our look back distance is finite 
to the visible universe, for any given I in the linear part we have a minimal q for 
which the integrand is still nonzero and below that we won't have any contribution 
to the power. This property shows up in the angular spectrum as a suppression of 
the power for the high Vs. for such Vs qmm is already in the region where the power 
is in the suppressed part of the transfer function. For the redshifts shown in Fig. |l| 
this suppression starts above / = 1000. It can be seen most clearly for the graph 
for z = 0.4. In fig. ^ we can see the effect for all three Vs at the very low redshifts. 
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10-^ I 

0.5 1.0 1.5 2.0 2.5 3.0 

z 

Figure 2: The shear spin-weight 2 angular power spectrum for different Vs as a function 
of z. The cosmological model is ACDM with {^de = 0.7, fift- = 0,n = \,wde = 
The three lines correspond to Multipols / = 10, 30, 300 from the bottom up, respectively. 

However, a range of red-shifts exists for which the spectrum's evolution is similar for 
all Vs and determined by $r as we will show next. 

4. red-shift evolution of the cosmic shear angular power spec- 
trum 

In this section we would like to describe and explain some qualitative aspects of the 
redshift evolution of 3D angular spectrum using our solution for the functional de- 
pendence of the shear on wtot{z). We will show that the 2Ci{zys depend weakly on 
the shape of the initial spectrum (the spectrum at the beginning of matter domina- 
tion) and more importantly that the redshift evolution is quite sensitive to changes 
in Wtot independently of the initial spectrum. Since we wish to focus on determining 
the qualitative aspects of the spectrum's redshift evolution, we will use approxima- 
tions that are simple enough to allow us to obtain analytical expressions. Of course, 
the spectrum can be evaluated numerically very accurately for each specific cosmol- 
ogy using CMBFAST/CAMB etc.. using the same techniques that are applied to 
the CMB. In addition to the qualitative analysis we will show accurate numerical 
solutions for the case of ACDM (using CAME) to supplement our qualitative results. 

In the above and subsequent calculations and analysis several assumptions are 
made for the sake of simplicity. The validity of these assumptions is generally ac- 
cepted. However, the sensitivity of our results to them should be explored further 
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either analytically or numerically. One such assumption that may be particularly 
relevant is spatial flatness. We hope to discuss it in future work. 

From the epoch of cold matter dominance onward, the time evolution of the 
perturbation $fc is independent of k (under the assumptions specifled previously). 
From the solution eqs.( |2.19| )-( |2.23| ) one can see that the time dependent part $t(-2) 
has a simple and known functional dependence on Wtot- It follows that if we had 
full knowledge of the 3D spectrum at different times we could extract Wtot- It was 
shown that the 3D spectrum in fc— space can be used to constrain the cosmological 
parameters |]52| , |28[| . We have found that Wtot (and thus the expansion history of 
the universe) can also be constrained from the redshift evolution of angular power 
spectrum multipoles. However, a subtle effect complicates matters: Looking farther 
in the radial direction necessarily involves looking back in time. This effect forces a 
mixing of the spatial and time dependent parts of the perturbation. 

The 3D spectrum P-yik) evolves in time independently of k. Its multipole ex- 
pansion in terms of the multipoles 2C1 is, however, distance-dependent and therefore 
the redshift evolution of the different multipoles becomes /c-dependent. We will show 
that this dependence on k is advantageous and useful. In the expression for the 
multipoles both $t(-2) and the luminosity distance dilz) appear. Both quantities 
depend on Wtot{z) in a different way making them more sensitive to changes in Wtot 
at different redshifts. More precisely, their sensitivity to changes in wtot varies in 
opposite ways. We have found that due to this, in the expression for the multipoles, 
their explicit combination is less degenerate in Wtot than each of them separately. 

4.1 The qualitative dependence of the spectrum on $7^ 

To understand how the 2Ci's depend on ^t{z) we will assume for the moment a 
flat (constant) initial spectrum. We discuss a flat spectrum because it gives us 
better insight to the qualitative behavior of the solutions. We will then relax this 
oversimpliflcation. A flat initial spectrum corresponds to assuming a flat primordial 
spectrum (n = 1) and Tq{zin) = 1. The solution in this case can be obtained 
analytically in a closed form. Moreover, with a constant transfer function in eqs. (|3.9[ - 
|3.1CI| ) it becomes possible to calculate the angular moments for any spectral index 
in terms of hypergeometric functions. For a flat spectrum (n = 1) the result is 
extremely simple 




(4.1) 



The function in eq. (|4.1|) is essentially a delta function. 




(4.2) 
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whose normalization is determined by the integral J^^ du2 {^j^ + '^'^2 ~ 
■jTZi- The approximation in eg. ( [4 .21) becomes better for larger /'s, exactly the range 
of /'s of interest. Putting everything together we obtain, 

2Q = (47r)M^|±|: [du, fdu, ""' - - ^t{u,)^t{u,)--^^u,5{u, - u,) 
{I -2)1 J J wiui W2U2 - 1) 



min[wi ,W2] 

= {A.fAl{l + 2) ! du {w,-u)iw2-u) 

J U W1W2 


We can assume without loss of generality that wi < W2, so our final result for this 
case is 

2Q{w„ W2) = {AnfAlil + 2) [ ^ 

J U W1W2 



We observe that the 2C1S depend on Wtot through an integral expression of the square 
of $T and a kernel function of m = . We recall that eg .(14.41) is derived assuming 



1+2 

a constant transfer function. Since $t(0) is finite, the integral in eg. ( |4.4|) is formally 
divergent. The formal divergence at small u is not physical rather it is a property 
of blue or fiat spectra (n > 1). In technical terms, tracing back the properties of 
the small u region of the integrand, one sees that it corresponds to the region of 
high g's in the integrand of eg.( |4.1| ). To further understand the behavior of our 
integral and the infiuence of the transfer function on the result let us look at the 
evolution of perturbations before the epoch of matter domination. During matter 
domination $ is constant but during radiation domination the solutions for $ inside 
the horizon decay as l/a? ~ l/??^- For values of q larger than some maximal value 
5'max the solutions of $ are therefore completely suppressed, and thus the u integral 
is effectively cutoff at m oc 1/gmax and becomes finite. 

The real transfer function undergoes a smooth transition from unity at small q 
to zero at high g's (it should be evaluated in the transition region by approximations 
such as BBKS). If we use accurate approximations of the transfer functions shape 
in Eg.(^]9D the analytic solutions become too complicated and it is hard to gain any 
understanding from it. Rather, we will demonstrate the effect of the suppression of 
high g's using a step function approximation of the transfer function, representing 
sharp cutoff on the spectrum at some q = gmax- 

When the transfer function can be approximated by a step function Tq = 1 — 
0{q — q) then the integral in eg. ( |2.23|) becomes 

00 g 

f —3i{uiq)ii{u2q){l -0{q-q})= [ — J?(gMi) ^ " "2)- (4.5) 

J q J q + 1 
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The definite integral 21(1 + 1) / ^ji'^{x) depends for large / on x/l so it can be 



approximated (numerically) for large /'s and for (x//) > 1 as 



where a,b are some /-independent constants. The integral vanishes for (x/l) < 1. 
The expression for the shear becomes 

,C, = 4A(2.n, + 2) / f ( Hl^) ( '^f^) <S.Au)%.7) 

J \ wi J \ W2 J \ ma J U 

The resulting expression ( [4.7D is different from the one in eq.( [4.4| ) in that its inte- 
grand has an added /-dependent numerical factor (|4.tj|) and, more importantly, it is 



finite. The functional dependence on wtot of both expressions comes from so this 
dependence is shared whether one assumes an unrealistic trivial transfer function 
{Tq = 1) or a more realistic step function form to the transfer function. Based on 
this comparison, we expect that the functional dependence of the solution with the 
correct, fully realistic and numerically calculated transfer function will also have the 
same simple dependence on Wtot that comes from 

4.2 The sensitivity of the 2^'^ evolution with red-shift to changes in Wtot 

In the previous subsection we have determined the dependence of the angular power 
spectrum on But, it is the combination of the evolution of the perturbations 
convoluted with the geometrical kernel that determines the true sensitivity of the 
angular spectrum to changes in Wtot- In Fig- i| we show the 2C1S for four different 
cosmological models. Also shown is the luminosity distance di-, and $r for each of 
the models. It is evident from the figure that the $t's are more sensitive to changes 
in the expansion history at low redshifts while the c/i(z)'s are more sensitive to such 
changes at higher redshifts. It is also evident that 2C1S sensitivity follows ^t, so 
we may argue that the 2C1S "measure" The sensitivity of the 2C1S is twice 
as much as that of $t simply because they depend on ($t)^- This means that to 
distinguish between these models using luminosity distance measurements we would 
need accuracies at the sub-percent level, while they could be easily distinguished if 
the 2C1S can be measured with about a percent accuracy. The realistic prospects of 
constraining the total equation of state and other cosmological functions depends on 
how well it will be possible to handle the low redshift range and on the selections that 
are made in the survey and in the data analysis. Hence, proper numerical analysis of 
the power spectrum's sensitivity ^ '^^^ and the likelihood calculations are necessary. 
These are left for future work. 
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Figure 3: The three panels show the relative difference (in percent) as func- 
tion of red-shift between different cosmological models and a ACDM concordance 
model. The four lines correspond to the four sets of {^de^wde) parameters: 
(0.72, -0.95), (0.71, -0.975), (0.69, -1.025), (0.68, -1.05). The first panel shows the rela- 
tive difference of the 2C1S for I = 100, the second is of $t and the third shows di. 
and differ in their sensitivity. 

5. Conclusions and outlook 

Several studies of the ability to use 3D weak lensing measurements to constrain cos- 
mological parameters exist in the literature. To the best of our knowledge, all these 
studies have presented the numerical results for specific fiducial lensing surveys with 
weakly coupled DE that has an equation of state of the form wde = Wq + Wiz/ {l + z) 
and positive speed of sound. These studies are important in demonstrating the ad- 
vantages of three dimensional analysis over the two dimensional one. However, it 
is not easy to extract from them information on the functional dependence of the 
3D spectrum on the equation of state and hence on the degeneracies with respect to 
simultaneous changes in several cosmological parameters. Hence, it is harder to use 
them to forecast the possible improvements in the constraints that one should get 
from more accurate shear measurements. 

By calculating the functional dependence of quantities such as the 3D cosmic 
shear on Wtot in a fully relativistic and without assuming a specific functional form for 
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Wtot, we can estimate the theoretical hmits on the constraints from a given accuracy 
for shear measurements. Our conclusion is that the red shift evolution of the spin 
weight two angular power spectrum of shear is the most sensitive measure. In this 
paper we have presented our analytical results and some numerical examples for the 
simple case of Linder's parameterization thus laying the ground for a more detailed, 
application oriented, analysis. We have found that redshift evolution of the angular 
power spectrum of the three dimensional cosmic shear can be used to determine the 
expansion history of the universe. This method is sensitive to other prior assumptions 
about the functional dependence of Wtot on redshift than kinematic methods that rely 
on the homogeneous and isotropic universe (such as luminosity distance etc.). 

The implementation of our method will require extraordinary observational ef- 
forts because it requires very accurate measurements of the shear. Further, good 
determination of the three dimensional cosmic shear requires accurate determina- 
tion of redshifts of the sources, and in addition a good measurement of luminosity 
distance (or angular distance) as a function of redshift. To achieve such ambitious 
observational goals would probably require combining the results of several accurate 
experiments. 

The equations for the metric perturbation $, unlike the equations for the density 
perturbation 6p, can be adapted in a straightforward manner to other metric theories 
of gravity which modify Einstein's general relativity. The equations can be easily 
modified to account for various corrections, and hence our methods can be used, in 
principle, to test modifications to general relativity. A discussion of this point can 



be found in 1153 
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A. Expressing Wtot in terms of $ 

We would like to transform eq.( |2.19| ) into a more convenient form. First, we change 
the integration variable from to p by using the conservation equation 

dlnp= -3{l + w)dlna. (A.l) 

Since dhia/drj = Tiijf), it follows that drj = d\na-^^-^ = — 1/3 ^^^^^^ , then the 
growing solution is given by 

<l>+(r/,f) = 4vrGC2(f)^ / \a'^ (A.2) 

a J 6 rip 

(the sign was absorbed in changing the integration boundaries.) 

If the universe is spatially flat then from Friedman's equation it follows that 
U'^ = ^a^p. In this case eq.([Al2D simplifies, 

*.fe2) = y^C.(f)^/<ip^. (A.3) 

Because the spatially fiat case is simpler we will assume a spatially fiat universe for 
the rest of this discussion. 
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We may further change variables to 

X = (A.4) 

Here Oq and po are the values of the scale factor and the energy density today. The 
value of the variable x today is Xo = 1 and it vanishes for very early times, if 
the universe was matter dominated as expected. Consequently, the range of x is 
0<X< 1. 

Using eq.( |A.l| ) we get d\nx = d\na — l/2d\np, so d\np = — ^^^i^^^^^^^filnx 
and therefore dp^^ = x = d\n p — 3(^1+^170+2 '^^ finally get 

<|.+(x,x) = V2A^C,{x)- jdx^^^P^y (A.5) 

XJ 5 + 3wtot(x) 



The solution x) factorizes 

=C(f)$T(x) (A.6) 



C{x) = V2AnGC2{x) (A.7) 



We can invert eq.(^), since (x$t) = j}^^^- It follows that 

l~5d^{x<l>T) . . 

If we so wish we can also express $r and Wtot as a function of redshift z = ao/a—l. 
Since x = x = corresponds to 2 ^ oo, and X = 1 corresponds to z = 0. For a 
spatially flat universe, as we are considering 

^ g/ao ^ 1 1 .A inx 



From eq. (|A.8|) we find 



'^t{z) = 4- fdz'^ '^^'^f^ , (A.11) 
Xiz) J dz 5 + 3wtot(i) ^ 



or equivalently, a differential equation for $t 
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From eq.([/OD it follows that d, \nx{z) = d, ln(a/ao)-i<9^ Inp, so = d, ln(Y^) + 



1 + 'Wtot)dz^^{j^), and the final result is that 



1 



31+ Wtotiz) 



X{z) l + z 2 1 + z 

Substituting eq.( |A.13| ) into eg. ( |A.12| ) we get 

2(1 + 2) 



5 + 3wtotiz 



-dMz) 



5 + Swtotjz 
2(1 + z) 



1 + ^'tot(^) 
5 + 3wtot{z 



(A.13) 



(A.14) 



Equation ( |A.14| ) allows us to solve for $r(z) in terms of Wtot- We may define an 
integration factor 

5 + Swtotiz) 



Kz) 



dz 



2(1 + z) 



(A.15) 



that simplifies eq.( |A.14j ) 



1 + Wtot (2:) 

2(1 + z) 



(A.16) 



The initial conditions on $t are prescribed at the initial redshift z^^. In terms of 

Kz), 



^t(z) 



Zin 

$T(..n) + e-^(^) jdz e^(^ (A.17) 



Here we have integrated eq. ( [A. 141 ) from large redshifts towards smaller ones. We can 
also integrate eq.( |A.14| ) from small redshifts towards larger ones. This will require 
knowing the amplitude of the perturbation at late times which is harder to determine. 
From eq. ([A.14|) we can also solve for Wtot in terms of $t, 



Wtot{z) 



2{1 + z)dz^T{z) + 1 - 5^t{z) 
1 - 3<1>t{z) 



(A.18) 



B. The relation between $r(^) and the hnear growth factor 

For the ACDM model there is an exact solution for the evolution of the matter 
density perturbations in the linear regime. The growth rate F{z)/(^1 + z) for the 



linear density perturbations 6m{z) = 6{0)F{z)/{l + 2:) is ||49 



oc (1 + z 



H{z) 



dz'(l + z') 



Ho 
H{z') 



(B.l) 
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To find the relation between F(z) and our solution for ^riz) we use eq.( 2.25 ) in the 
Newtonian approximation 



5p 



so that 



p 3{aHy 

Spm 5p p 



Pv 



p pm 3Ho^ ^moil + Z)' 



(B.2) 



(B.3) 



Thus F{z) OC <^T{z)/nmo. 

To show that the solutions are indeed identical we start from the solution for 
the metric perturbation. The growing mode solution of eq. ( |2.19D in cosmic time is 
given by 



$T(t) = ^ / diaii){l+wtotii)). 



(B.4) 



Using the fact that = ^^p and Hdt = —dz/{l + z) we express $t as a function 
of z 



^t{z)^{1 + z)H{z) J — 



-^^^——{1 + Wtotiz')) 

zYH[z') 



[l + z)H{z) 



dz' 1 + z' 
H\z') (1 + ^')^ 



(B.5) 



ForACDM, wde = —1 and therefore p + p = Pm- Using the fact that p^ = Pm(0)(l + 
z)^ we finally get 



<!>t(z) OC Qrnnil + Z 



H{z) 
Ho 



dz'il + z') 



Ho 



H(z' 



(B.6) 



where = ^m(O)- The solution in eq. ( [B.ti| ) is identical to the one appearing in 
eq. (IB.ll) . The growth factor in eq. (|B.1|) is normalized such that -F(O) = 1. The 
normalization that we choose in the paper following is such that ^t{z ^ 1) = 1/5 
outside the horizon and thus numerically 



mo 



T 



l + z) 



Mz) 

Ho 



dz'il + z') 



Ho 
H{z') 



(B.7) 



The full amplitude of the metric perturbation is related to the density perturbations 
amplitude through the same relation as in ||^ . 
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